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Abstract. Photonic quantum simulators are promising candidates for providing 
insight into other small- to medium-sized quantum systems. The available photonic 
quantum technology is reaching the state where significant advantages arise for the 
quantum simulation of interesting questions in Heisenberg spin systems. Here we 
experimentally simulate such spin systems with single photons and linear optics. The 
effective Heisenberg-type interactions among individual single photons are realized by 
quantum interference at the tunable direction coupler followed by the measurement 
process. The effective interactions are characterized by comparing the entanglement 
dynamics using pairwise concurrence of a four-photon quantum system. We further 
show that photonic quantum simulations of generalized Heisenberg interactions on 
a four-site square lattice and a six-site checkerboard lattice are in reach of current 
technology. 
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1. Introduction 

More than a quarter of a century ago, Richard Feynman [l]|2] envisioned that a 
weU-controUed quantum mechanical system can be used for the efficient simulation 
of other quantum systems and thus would be capable of calculating properties that 
are unfeasible for classical computers. Quantum simulation promises potential returns 
in understanding detailed quantum phenomenon of inaccessible quantum systems, 
from molecular structure to the behavior of high-temperature superconductors |3,4|. 
Moreover, quantum simulations are conjectured to be less demanding than quantum 
computations by being less stringent on explicit gate operations or error correction fi] . 

Due to these reasons, quantum simulation has led to many theoretical proposals |4j- 
[s]. Recently, various quantum simulators based on different physical platforms are 
being constructed, such as atoms in optical lattices [9-14 , trapped ions 15-20 , nuclear 



magnetic resonance systems [21|[22|, superconducting circuits [23 1, as well as single 
photons 24-32 . Motivated by the seminal work of Knill, Laflamme and Milburn |33| 
photons are proven to be a suitable system for efficient quantum computing and quantum 
simulation. Precise single-photon manipulations and tunable measurement-induced 
two-photon interactions are the essential ingredients for photonic analog quantum 
simulation and have been demonstrated. In addition to high-level quantum control, such 
photonic quantum simulators can produce exotic entangled states which are important 
for understanding the many-body dynamics in quantum chemistry and solid-state 
physics [7(|8m. 

Photons barely interact with the environment, which makes them ideal for 
exploiting various quantum phenomena. But in contrast to ions or atoms in optical 
lattices where physical interactions can be easily implemented via Coulomb interactions 
or Feshbach resonances, photonic simulations require various preparation steps to 
achieve controllable effective interactions. First, instead of non-correlated particles, 
one creates entangled singlet pairs of photons from non-linear crystals and these are 
used as input states for quantum simulation. In the language of quantum information 
we would say that each photon brings half an ebit (one half of the singlet) of correlation. 
Second, the non-correlated photons are overlapped at a tuneable directional coupler and 
measured in different output modes to introduce effective non-linearities that create the 
targeted many-body correlation. This is an effective interaction among the photons (for 
detailed calculation see Ref. |30|). Note that the recently proposed universal quantum 
computation based on coherent photon conversion provides an alternative avenue for 
achieving photon-photon interaction |34|. 



2. Analog quantum simulations with photons and linear optics 

In the case of the simulation of spin- 1/2 particles the photon's polarization is 
ideally suited as horizontally-polarized states \H) and vertically-polarized states \V) 
representing for example, spin-up and spin-down states. Moreover, the ability to 
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prepare symmetric polarization-entangled states |$)^2 = -^{\H)i\H)2 ± |l^)i|^)2), 
I*)t2 = 73(1^)11^)2 + 1^)11^)2) and the anti-symmetric state |^-)i,2 = ^{\H)i\V)2- 
\V)i\H)2) enables the establishment of states with bosonic and fermionic character [29[ 
30,35] . The latter shares the same quantum correlations as Heisenberg-inter acting spins 
or so-called valence bond states [36j . 

The theoretical investigation of strongly-correlated spin systems has led to few exact 
theorems which in some cases are of importance for the quantum simulation of chemical 
and physical models. In the particular case of a nearest-neighbor antiferromagnetic 
Heisenberg-interacting spin system it was shown by Marshall [37] that the ground 
state for N spins on a bipartite lattice has total spin zero (aS^ = 0). This theorem 
and its extension |38| lead to the fact that the ground state must be built as a linear 
superposition of singlet spin states or valence bonds. This constraint, that forces the 
ground state's total spin to be zero, gives rise to various valence-bond configurations 
that are either localized or fiuctuating as superposition of different singlet partitionings. 
Localized configurations are typically referred to valence bond solids and delocalized 
valence-bond states correspond to frustrated quantum spin liquids or resonating valence- 



bond states 39,40i|. Recently, the photonic quantum simulation of a four spin- 1/2 square 



lattice as an archetype system [30] showed that quantum monogamy [41|[42] plays an 
important role in frustrated Heisenberg spin systems. 

Analog quantum simulation uses the advantages of digital and adiabatic quantum 
simulators by utilizing continuously tunable quantum gates. In Ref. |30|, we reported the 
analog quantum simulation of the ground state to probe Heisenberg-type interactions of 
a spin-1/2 tetramer. The variable measurement-induced interactions between photons 
are crucial for our experiments. In conjunction with single-photon manipulation and 
detection, it allows us to show the dynamics of ground-state energies and pair-wise 
quantum correlations of this tetramer. Here, we present a detailed experimental study 
of these photon-photon interactions. Depending on the interaction strength, frustration 
within the system emerges such that the ground state evolves from a localized to a 
delocalized valence-bond state. The quantum monogamy relation is unambiguously 
demonstrated via entanglement measures [4T, 



42 



The experimental setup for studying variable measurement-induced interactions is 
shown in Fig. [TJ Our pump source is a mode-locked Ti:sapphire femto-second laser with 
a pulse duration of 140 fs and a repetition rate of 80 MHz. The central wavelength of 
the pump is at 808 nm. Then we use a /3-barium borate crystal (BBOO) to up-convert 
the pump pulses to ultraviolet (UV) pulses via second harmonic generation. The up- 
converted UV pulses' central wavelength is 404 nm with a power of 700 mW. Then we 
clean the UV pulses with several dichroic mirrors (DM). Photons 1 and 2 are generated 
from a BBO crystal (BBOl) via spontaneous parametric down conversion (SPDC) in 
a non-collinear type-II phase matching configuration and in an polarization-entangled 
state after walk-oflF compensation |43|. Photons 3 and 4 are generated from BB02 in a 
collinear type-II phase matching configuration and are separated by a polarizing beam 
splitter (PBS). We guide photons 1 and 3 to a tunable directional coupler (TDC) and 
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Figure 1. Experimental setup for testing variable measurement-induced interaction 
between two independent photons. Ultraviolet (UV) femtosecond laser pulses are up- 
converted from a mode-locked Tiisapphire laser. The laser pulses are up converted 
(BBOO) and the obtained temporal duration and repetition rate is 140 fs and 80 MHz, 
respectively. The UV pulses and the remaining fundamental pulses are separated with 
several dichroic mirrors (DM). Lenses are used to focus and collimate the UV pulses in 
order to achieve good up and down conversion efficiencies. One pair of the polarization 
entangled photons (photons 1 and 2) is generated from a /3-barium borate crystal 
(BBOl) via spontaneous parametric down conversion (SPDC) in a non-collinear type- 
II phase matching configuration. Two half-wave plates (HWP) and compensating BBO 
crystals are used to counter walk-off effects in the down-conversion crystal. Another 
pair of correlated photons (photons 3 and 4) is generated from BB02 in a collinear type- 
II phase matching configuration and are separated by a polarizing beam splitter (PBS). 
Single-mode fibers and interference filters (IF) are used to clean their spatial and 
spectral modes, respectively. We vary the path length difference between two photons 
with a motorized translation stage mounted on the fiber coupling stage of photon 
1. Fiber polarization controllers (PC) are employed to eliminate the polarization 
distinguishability of the two interfering photons. The tunable measurement-induced 
interaction among photons is achieved by a tunable directional coupler (TDC), followed 
by a projective measurement of one photon in each of the four output modes. 



then detect them by avalanche photodiodes (APD), which together enable the tunable 
measurement-induced photon-photon interactions to happen. The relative temporal 
delay between the photons is adjusted with a motorized translating stage mounted on 
the fiber coupler of photon 1. Fiber polarization controllers (PC) are employed to 
eliminate the polarization distinguishability of the two interfering photons. 

Bunching due to the two-photon Hong-Ou-Mandel (HOM) interference |44| and 
the corresponding anti-bunching effect |35| are crucial for many quantum information 
processing protocols, especially for photonic quantum computation experiments (C- 
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phase gate [45}|46], entanglement swapping 47,48 , etc.), as well as for our photonic 
quantum simulation f30]. The bosonic nature of the photons shows up when the input 
states are superimposed at the TDC such that a detection, even in principle, cannot 
distinguish either of them. This leads to a superposition of double occupations on both 
outgoing modes and thus suppression of the coincidence detections, where one photon is 
detected in each output mode. The visibility of this HOM dip is one when the TDC is 
set to have equal splitting of transmitted and reflected photons similar to a 50/50 beam 
splitter. As soon as the two input photons can be partially distinguished by unbalancing 
the splitting ratio the visibility decreases. The dependence of the ideal visibility {Videai) 
upon the reflectivity of TDC {rf) is the following: 

1 - 2ry + 27y2 

In Fig. [2] this reflectivity dependent ideal visibility is plotted in the black solid curve. 
Experimental imperfections due to high-order emissions from SPDC and group velocity 
mismatch reduce the visibility. Fig. [2] also shows the measured visibilities (black squares) 
and their corresponding fit (dashed red curve). 

3. Pair- wise entanglement dynamics 

The main advantage of the precise quantum control of individual particles is that inter- 
particle entanglement dynamics can be investigated. By using a similar experimental 
configuration as in Ref. |30|, we study the entanglement distributions among different 
particles with respect to the effective interaction strength that was tuned by the TDC. 
For the quantification of the bipartite entanglement in our system, we use the measure of 
concurrence [49], which, for a given state p, is C{p) = max{0, ^/)^ — ^/)v2 — ^/\^ — ^/\^} ^ 
where are the eigenvalues of the matrix pTip^Ti in non-increasing order by magnitude 
with S = cTy CTy, where ay = —i\0){l\ + i|l)(0|. While the previous quantum 
simulation characterized the pair- wise energy dynamics of Heisenberg interactions that 
were directly extract able from measured coincidence counts, this experiment requires 
the reconstruction of the density matrices to obtain the concurrence values. In the 
experiment we tune the refiectivity of the TDC and hence vary the photon-photon 
interaction strength. Various four-photon quantum states are tomographically measured 



and the density matrices of them are reconstructed 50,51 . The concurrence of the two- 
photon subsystems is calculated from the four-photon density matrices by tracing out 
the other two photons. 

Due to the quantum monogamy relations (30[[4l}|42] the total amount of pair- 
wise entanglement stays constant while the change of interaction strength affects the 
distribution and thus the ground state configurations. When tracking the change in 
entanglement by using concurrence, the expected "sudden death" and "sudden birth" 
of entanglement l53i can be seen (Fig. Isl). While this concept is typically used for 



studying environment-induced decoherence 54 , similar behavior can be observed here 



too. In fact our tunable interactions allow to mimic a controlled interaction with the 
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Figure 2. Experimental demonstration of variable measurement-induced interaction 
between two independent photons, (a), (b) and (c) Measured four- fold coincidence 
counts are plotted versus the relative optical delay between the interfering photons 
and fitted with a Gaussian function, when the reflectivities of the TDC are tuned 
to be 0.17, 0.5 and 0.67, respectively, (d) Visibility of the Hong-Ou-Mandel (HOM) 
dip when measuring four-fold coincidences using a TDC. The black squares are the 
experimental results. The black solid curve is a theoretical prediction based on Eq. [l] 
and the red dashed curve is a fit with the only free parameter of Vsys^ which is about 
0.853. The main reason of this non-ideal value is due to the higher-order emission 
generated from SPDC sources. The error bars are based on a Poissonian statistics. 

environment of two additional particles, which opens the possibility to obtain insights 
into complex decoherence mechanisms. We explicitly show the concurrence of different 
photon-pair configuations with respect to the TDC angle 9. The relation of the TDC 
angle and its reflectivity is given by 9 = arctan y^. One can see that the concurrence 
for one photon pair, e.g. Cu decreases rapidly as we increase 9. Around 9 = 0.274, 
all of a sudden the entanglement between another spin pair (e.g. Cis) is born at the 
cost of the reduced concurrence Cu^ which vanishes as 9 is further increased. The 
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Figure 3. Measured concurrence for various ground state configurations. 

Entanglement distribution between photons 1 & 2 (black square, fitted with black 
curve), 1 & 3 (red circle, fitted with red curve) and 1 & 4 (blue triangle, fitted with blue 
curve) are shown as a function of the TDC parameter 0. The rise of entanglement in 
one pair causes the fall of the entanglement in other pairs according to monogamy. The 
monogamy gives rise of the entanglement sudden death and sudden birth 52- 54 . We 



derive the uncertainties in concurrence from the density matrices, which are calculated 
using a Monte Carlo routine and assumed Poissonian errors. 



observation of similar disappearance or emergence of entanglement among the other 
photons demonstrates the capabihty of our quantum simulator for manipulating the 
quantum correlations (52] - [54| . 

4. Generalized Heisenberg spin model on a four-site square lattice and a 
six-site checkerboard lattice 

The fact that the properties of photons seem to make the simulation of ground 
states for Heisenberg-interacting spins feasible, opens promising perspectives for the 
investigation of more complex interactions by using current technology. Therefore we 
investigate generalized Heisenberg models for the four-site square lattice and the six-site 
checkerboard lattice. For the four-site square lattice we extend the model in Ref. |30| 
by adding a next-nearest neighbor interaction term (Fig 4a). The Hamiltonian for this 
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Figure 4. a) A four-site square lattice with Ji being horizontal coupling, J2 being 
vertical coupling and J3 being diagonal coupling, b) A six-site checkerboard lattice. 
Ji is the coupling between the nearest neighbors and J2 is the coupling on the cross 
bonds. 
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Figure 5. Ground state phase diagram for various parameter regimes. The 
coefficients of the ground state, \^g) = a\^^) + /^|*^*||), are shown for different ratios 
of J2/J1 and J3/J1. In a), b) and c), we show \P\ and |a| + respectively. A 
clear phase transition can be seen along the dashed line a) for Ji = J3, b) for J2 = J3, 
and c) for Ji = J2. 



system is: 

H = Ji(Si • S2 + S3 • S4) + J2(Si • S3 + S2 • S4) + J3(Si • S4 + S2 • S3) 

We consider the antiferromagnetic case with the couplings Ji, J2, Js ^ 0. For 
each of the three terms in i/, the ground state is a pair of singlets, \^=) = 1^^12)1^34)5 



1*11) = 1*13)1*24) and |$x) = 1*14)1*23) = l*=) - l*||)- In Ref- 30 we showed that by 



tuning J2/-^i from to oc, the ground state, gradually changes from \^=) to |$||). 

With the introduction of J3, the ground state of the system is still a superposition 
1$^) = a\^=) + with normalization condition 2{\a\^ + + \a + /3\^) = 1. 

Remarkably, tuning J3/J1 can induce sharp phase transitions with a sudden 
change of the ground state configuration due to the competing of the valence bond 
configurations. In Fig. 5, the ground state configurations for different coupling regimes 
are shown. There are three particular interesting phase transitions: 

• For Ji = J3, 1$^) suddenly changes from \^=) + |$x) to |$||) when J2/-^i is tuned 
across 1 
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Figure 6. Energy spectrum of the low lying eigenstates of the six-site checkerboard 
lattice system. The energy levels are shown as a function of the ratio J2/J1. For the 
ratio J2/J1 = 1 an avoid level crossing of the ground state occurs. 

• For J2 = J3, \^g) suddenly changes from \^=) to |$||) — |$x) when J2/J1 is tuned 
across 1 

• For Ji = J2, \^g) suddenly changes from + |$||) to |$x) when J3/J1 is tuned 
across 1 

The case with Ji = J2 7^ Js is widely studied for square lattice systems due to its 
relevance to cuprates, Fe-based superconductors, and other materials [55|[56|. Previous 
studies have shown that in the thermodynamic limit, when J3 > Ji, the system is 
in a diagonal Neel ordered state [58], which is consistent with the ground state |$x) 
for a minimum of four sites, as discussed above. In the regime where J2 < Ji the 
configuration for the ground state is still under debate due to regions that appear to 
be non magnetic. Numerical calculation have recently shown that this region is highly 
likely to be a quantum spin liquid with Z2 topological order [58] . However, for only four 
spins, the region with J2 < Ji has only a single ground state configuration + |$||). 

We also investigate theoretically the ground states of a Ji — J2 Heisenberg model on 
a six-site checkerboard lattice. The geometry of this system is shown in Fig. [4|3, where 
Ji is the coupling strength between the nearest neighbor sites and J2 is the coupling 
strength on the cross bonds. The coupling ratio J2/J1 is the only parameter of this 
system. The introduction of the next nearest neighbor coupling J2 makes this system a 
simple frustrated magnetic model. The study of this model on thermodynamics limits is 
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Figure 7. The independent ground state configurations of the six-site checkerboard 
lattice. The four different dimer coverings are labeled as a)|?/^i), h)\ip2)^ c)|?/^3), and 
d)|?/^4) in our system. 



motivated by the three-dimensional pyrochlore materials 59,60 . The two-dimensional 



model has been studied by several groups (61-67 . It is known that for the regime where 



the coupling ratio J2/J1 <C 1 the system has a co-linear Neel order. At J2/J1 ^ I 
numerical calculation suggests a plaquette valence bond solid ground state while 
the ground state at J2/Ji^ I remains under debate j^^,^. For this reason future 
photonic quantum simulations might provide answers to these open questions. 

Numerical studies of the six-site checkerboard system were done by diagonalizing 
the Hamiltonian with open boundary conditions. In Fig.|6]the energy spectrum of the six 
low lying states as a function of the ratio J2/J1 in the = subspace is presented. 

The energy spectrum shows an avoid level crossing around J2/J1 ^ 1, which indicates 
a dramatic change of the ground state properties there. In analogy to the Ji — J2 
Heisenberg model on a plaquette, where the ground states can be expressed as linear 
superposition of two different dimer coverings whose coefficients depend on the ratio 
J2/Ji-> the ground states for the six-site checkerboard lattice can also be expressed as 
superpositions of various dimer configurations. For the discussed six-site lattice system, 
fifteen different dimer coverings exist. However, only six are independent. By taking 
the symmetry of our system into consideration only four out of the six coverings are 
allowed (Fig.|7|). 

Therefore, the ground states can be described as superpositions of four dimer 
coverings, l^/^i), |'02), iV^a), |'04) in all region of J2/ Ji- In Fig. [s] we show the contribution 
of each dimer configuration with respect to the coupling ratio J2/ Ji- At J2/ Ji = 1 the 
coefficients for and 1^03) are equal and coefficients for \ip2) and 1^04) are exact zero. In 
our current convention of the dimer covering wave function this particular superposition 
gives us a plaquette state on the right four sites as it is shown in the inset. This coincides 
with the plaquette valence bond solid state in an infinite system. Remarkably, the six- 
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Figure 8. The six-site checkerboard lattice ground state. The according 
superposition of the four independent dimer configurations, IV^i), \i^2)^ I '^3)7 l'^4) is 
shown as function of the ratio J2/J1' The inset shows the plaquette ground state at 
J2/J1 = 1. 



site checkerboard lattice thus provides already a valuable hint for the true ground state 
in the thermodynamic limit. The ground states at the ratio J2 > Ji have close to equal 
contribution from \ipi) and 1^02), which suggests a significant contribution from a cross- 
dimer state on the left plaquette. We would like to mention that a cross-dimer ground 



state has also been suggested as a potential ground state for large J2 couplings 63 . It 



is interesting that only a small contribution from the dimer configuration {ip^} can be 
found in all the possible cases. 



5. Summary and Outlook 

In conclusion, today's available photonic quantum technology is reaching the stage 
where significant advantages arise for the simulation of particular interesting questions 
in solid-state physics and quantum chemistry. Therefore photonic quantum simulations 
provide exciting opportunities to cover, for example, the direct construction of custom- 
tailored many-body wave functions. Impressively, the usage of optical elements such 
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as tunable non-polarizing or polarizing beam splitters enables entangled few-photon 
states to construct many-body valence bond wave functions for molecular and solid- 
state systems due to non-classical interferences. As we have shown above, this is of 
particular interest in condensed matter physics as it provides insight into the frustration 
of strongly-correlated spin systems and the onset of quantum phase transitions. On the 
other hand, in quantum chemistry it effectively allows studying delocalized bonds in 
chemical structures and chemical reactions |5|. Thus being able to monitor the full 
dynamics of individual particles and bonds provide some fascinating perspective for the 
quantum simulation of small molecules or reactive centers. 

The main future challenge will be to increase the number of photons or degrees 
of freedom to realize a sufficient amount of qubits such that quantum computers can 
outperform their classical counter parts. In general, up to the level of approximately 
twenty qubits it presently appears possible to conceive a system based on bare physical 
qubits. However, given the current experimental limitations, operational fidelities and 
noise sources, it seems that useful system consisting of more than twenty qubits could not 
be realized without some level of error correction. But in contrast to the implementations 
of well-known quantum algorithms, such as Shor's algorithm for a computationally 
relevant key-length, the requirements for fault-tolerance are much less demanding. How 
much error correction is needed to achieve a useful quantum simulation is an open 
problem. 

Recent work |5y26j has shown that quantum systems with less than a dozen physical 
qubits are capable of simulating chemical systems with a precision that cannot be 
achieved by conventional computers, when processed via almost a thousand discrete 
gate operations. Although such small quantum systems are feasible by using present 
quantum technology, the requirements in terms of gate operations is tremendous. Thus, 
analog quantum gate operations look promising in reducing the technical complexity of 
performing such quantum simulation experiments by requiring fewer number of physical 
gates. 
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